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Abstract 

Neutralino proton cross sections are examined for models with R-parity 
invariance with universal soft breaking (mSUGRA) models, nonuniversal 
SUGRA models, and D-brane models. The region of parameter space where 
current dark matter detectors are sensitive, i.e. 1 x 10 -6 pb, is examined. 
For mSUGRA models, detectors are sampling parts of the parametr space for 
tan/3 ~ 25. The nonuniversal models can achieve cross sections that are a fac- 
tor of 10-100 bigger or smaller then the universal one and in the former case 
sample regions tan/3 ~ 4. The D-brane models considered require tan/3 ~ 15. 
The inclusion of CP violating phases reduces the cross section by a factor 
of ~ 2-3 (but also requires considerable fine tuning at the GUT scale). The 
expected particle spectra at accelerators are examined and seen to differ for 
each model. Three new regions of possible coannihilation are noted. 

I. INTRODUCTION 

Supersymmetric theories with R-parity invariance predict the existence of a dark mat- 
ter candidate, the lightest supersymmetric particle (the LSP). Accelerator and cosmological 
constraints then imply that the LSP is the lightest neutralino, x®, and that this particle 
is mainly gaugino ||. For heavy nuclei, x'j'-nucleus scattering is dominated by the spin 
independent part of the amplitude where the neutron (n) and proton (p) amplitudes are 
approximately equal. This allows one to extract from the data the spin independent cross 
section cr^o_ p . Current dark matter experiments are now sensitive enough to probe a signif- 
icant part of the SUSY parameter space. Thus DAMA and CDMS are sensitive to cr^o_ p in 
the range of (1 — 10) x 1CU 6 pb, and there will be perhaps a factor of 10 improvement in 
the near future. We ask here then what part of the SUSY parameter space can be tested 
for the range 

0.1 x 10" 6 pb < a^_ p < 10 x 10~ 6 pb. (1) 



We do this by examining the maximum theoretical cross section that lies in this domain 
by varying the SUSY parameters (e.g. tan/3, rn^o, etc). In the following, we consider three 
supergravity models: mSUGRA with universal soft breaking parameters [|3|-4Tl"|j, SUGRA 



models with non-universal soft breaking in both the Higgs 
squark/slepton sector 
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12|-|T7H and third generation 
and a D-brane model which allows for a specific pattern of 



nonuniversal gaugino, squark and slepton masses [|T^ . 

In our analysis ||19|| , we let the SUSY soft breaking parameters m g (gluino mass), m 
(scalar mass), A Q (cubic soft breaking mass) and tan/? =< H 2 > / < Hi > range over the 
following domain: 



rrig, n^o < 1 TeV; |y4 /m | < 5; 2 < tan(3 < 50. 



(2) 



In order that the analysis be accurate for this domain, we include the folllowing in the 
calculations; (i) we run the full 1-loop renormalization group equations (RGE) from the GUT 
scale M G = 2 x 10 16 GeV to the t -quark mass m t =175 GeV, iterating to get a consistent 
SUSY mass spectrum for a fixed set of GUT scale mass parameters, (ii) L-R mixing in 
sfermion mass matrices are included (which is important for large tan/3 for the third squark 
and slepton generations). (iii)One loop corrections to the Higgs mass matrix is included, 
(iv) One loop corrections to the mass rrib is included (which is also important for large tan/5 
where it produces a significant correction to the relation between (the b- Yukawa coupling 
constatnt) and m b . (v) QCD RGE corrections are included for contributions dominated by 
light quark masses, (vi) leading order (LO) and approximate NLO corrections [P0"1,P1 to the 
b — > S7 decay rate are included. We do not impose b — t unification at Mq (as is done in 
|16|), as this constraint is sensitive to possible unknown GUT scale physics. 

Various LEP, Tevatron and CLEO accelerator bounds limit the SUSY parameter space. 
For the light chargino Xi an d light higgs h we require 



m-± > 94GeV; m h > 95GeV, 

Al 

and for the B — > X s 7 branching ratio we use [^1 ] : 

1.8 x 1(T 4 < B(B -> X sl ) < 4.5 x 1(T 4 . 



(3) 



(4) 



The b-quark mass used is mb(mb) =4. 0-4. 5 GeV [22| and Tevatron bounds |23[] on m~ g and 
m,q are imposed. 

The theoretical analysis determines the Xi ~ Q scattering cross section, and in order to 



relate this to the proton cross section two parameters a and a^N enter [^J,^], where 

1 



[m u + m d )(p\uu + dd\p) 



(5) 



cr = ~{m u + m d )(p\uu + dd- 2ss\p) . 



65 MeV, in accord 
as 



as well as the quark mass ratio r = m s / \{m u + m^). We use here ovjv - - ' 
with recent analysis [E6| and a® = 30 MeV |25|]. The quark mass ratio is given in pT 
r=24.4±1.5. 

Theory allows one to calculate the mean Xi re hc cold matter (CDM) of the universe, 
parametrized by Q^oh 2 , where Q^o = p^o/p c (Px? = mass density of x?, Pc = 3Hq/8ttGn, 
H = Hubble constant) and h=H / (100km/ sMpc)). The combined world average for H 
is [^8| H Q = (71 ± 3 ± 7) km s _1 Mpc -1 . A number of astronomical analyses of matter 
density fl m gives |29"1,|30| Q m ~ 0.3, and using the baryon density Qb = 0.05, we assume that 
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fi^o = 0.25 ± 0.10 . In view of possible systematic errors that may exist, we chose here a 2 

0.02 < O^o/i 2 < 0.25 . (6) 



std range for Q^oh 2 : 



The theoretical formula for 1S gi ven by [31]: 

/T-o\ 3 / T \ 3 /V 1/2 

^o/i 2 = 2.48 x 10- 11 (-M (-^-) Xf / (7) 

Xl \Tj \2.73j dx(a ann v rel ) 1 ' 

where T^o is the freeze out temperature Tf{xf = kTf/m^o), Nf is the number of degrees of 
freedom at freezeout, (T^o/T 7 ) 3 is the reheating factor, a ann is the annihilation cross section, 
v re i is relative velocity, and ( ...) means thermal average. 

Eq.(6) puts a strong constraint on the SUSY parameter space and this then affects the 
Xi~ p cross section. From Fig. |l]one sees that (i) <J ann decreases with increasing mo and m^o 
(thus increasing Q^oh 2 ) and (ii) if 2m^o is near (but below) rrih, ma, mA rapid annihialtion 
occurs through s -channel pole (thus decreasing Q^oh 2 ). While LEP has now eliminated 
most of the parameter space where 2m^o ~ m^, both H and A become light for large tan/9 
making this effect latter significant in that domain. 

In general one sees that the bounds of Eq. (6) exert a strong inflence on the allowed 
SUSY parameter space, and this in turn affects the size of the Xi ~ 1 cross section. From 
Fig. |2| one has that (i)cr^o_ p becomes large for light (first generation) squarks (i.e. small m ) 
and light Higgs bosons, which is just the region of rapid early universe annihilation. Thus 
the lower bound on Q^ah 2 can produce an upper bound on cr^a_ p . (ii) We also note that 
<r^o_p increases with tan/3. 



II. MSUGRA MODELS 

The mSUGRA model with radiative breaking of SU(2) x U(l) depends upon four pa- 
rameters and one sign. These may be taken as following: mo, the universal scalar mass at 
Mg; mi/2, the universal gaugino mass at Mq (or alternately one may use m^o or nij at the 
electroweak scale since these scale with m.\/2 i.e. m^o ~ 0.4mi/2 and m~ g ~ 2.8m!/ 2 ); Ao, the 
universal soft breaking mass at M G ; tan/3 =< H 2 > / < Hi >; and the sign of /i, the Higgs 
mixing parameter in the superpotential (W^ = fiHiH 2 ). (With our choice of sign in W^, 
the b — > 57 constraint eliminates most of the /i > parameter space.) 

We proceed by varying the above parameters over the allowed space subject to all the 
constraints discussed above. The maximum value of cr^ -p as a function of m^o is given in 
Fig. |3|. We see that current experiments sensitive to o"^o_ p > 1 x 10 -6 pb are probing only 

large tan/3 regime for this model, i.e. tan/3 ~ 25. 

To obtain further insight as to the parameter space being probed by current experiments, 
we show in Fig. |] the value of f^o/i 2 when cr^o_ p takes on its maximum value for the 
characteristic case of tan/3=30. We see that f^o/2 2 increases as m^o increases (as expected 
from the general discussion above) from roughly the minimum to maximum value allowed 
by Eq.(6). Thus an accurate determination of Q m h 2 as might be expected from the MAP 
sattelite, would significantly help narrow the allowed SUSY parameter space. 
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One notes in Fig. 3 that the larger tan/3 cross section sustain with increasing m^o more 
than the smaller tan/5 do. This is due to the fact noted above that m# and uu become 
lighter as tan/3 increases increasing cr^o_ p . This is shown in Fig. |5] for the two cases of 
tan/3=30 and tan/3=50. In contrast, m.^ is relatively heavy, i.e. for tan/3=30 one has that 
m/j increases monotonically with m^o from 112 GeV to 128 GeV. Since mo is not large when 
&y-°-p i s a ^ maximum, the squarks lie below and close to the gluino where m- g ~ 7m^o. 



III. NONUNIVERSAL MODELS 



The possibility of nonuniversal soft breaking significantly changes the region of parameter 
space being accessed by current dark matter experiments. We consider here the case where 
nonuniversal soft breaking masses are allowed both in the Higgs and third generation sectors. 
A general parametrization at Mq then 



m k = m o(! + S h 
m q L = m o(! + h) 
m bl = m o(! + K 



rn 



2 _ _2 

h 2 



ml(l + 5 2 )- (8) 



m u R =m (l + 8 4 ); m eR = m (l + 5, 
m il = m 2 (l + 5 7 ). 



Here m is the universal mass of the first two generations, and 5i are the deviations for the 
Higgs and third generation. We assume here that — 1 < 5i < 1. (Note that for SU(5) one 
would have £ 3 = 64 = £5 and 5§ = 67.) 

In order to understand the effects that occur in this more complicated situation we first 
note thate Xi is a mixture of gaugino and higgsino part : Xi = ocW$ + (3B + 7^1 + 5H 2 . 
Further, the spin independent cross section arises due to the interference between the gaugino 
and higgsino parts of Xi> causing o"^o_ p to increase with increasing interference. The SUSY 
parameter that to a large extent controls the amount of interference is fi 2 , interference 
increasing (and hence cr^o_ p increasing) as /i 2 decreases, and interference decreases as /i 2 
increases. The value of fi 2 at the electroweak scale in terms of the GUT scale parameter is 
determined by the RGEs. In general, one must solve these numerically. However, one can 
get a qualitative understanding of the effects of the 5i from an analytic solution which is 
valid for low and intermediate tan/3 |14 |: 



V = ^23j[( 2 t 2 ^ (~2 — ^ 3 + 6i) — 6 2 + ^)\ m o (9) 

+ universal parts + loop corrections. 

Here t = tan(3, and D ~ 1 — ( y m t /2QQsinf3) 2 . In general D is small i.e D < 0.23. Eq. (9) 
shows that it is necessary to consider both the squark nonuniversalities (£3 and #4) as well 
as the Higgs (5i and 5 2 ) since they produce effects of comparable size. 

One sees from Eq. (9) that fi 2 will be significantly reduced (and hence cr^o_ p increased) 
if one chooses <5 3 , <5 4 , 5\ < and <5 2 >0. (The reverse will be the case for the opposite 
choice of signs.) This effect can be seen in Fig. ^, where the maximum value of o" x o_ p is 
plotted for tan/3 = 7 for the nonuniversal model (upper curve) and mSUGRA (lower curve). 
One sees that with the choice £3, £4, 5\ < and 5 2 >0 one can increase the cross section 
by a factor of 10 to 100. As a consequence, there are regions of parameter space where 
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the nonuniversal models can be probed to much lower tan/?. This is exhibited in Fig. [7]. 
where the maximum cross section is given for tan/?=7 and tan/?=5. One sees that current 
experiments are probing tan/3 as low as tan/3 ~4. Since these experiments have also excluded 
the region cr^o_ p ~ 10 x 1CT 6 pb, one finds that part of the parameter space for tan~ 12 
is already experimentally excluded for these models. Of course, the reverse choice of signs 
for Si will lower the cross section for any value of tan/?, leaving such parts of the parameter 
space mostly as yet unexplored experimentally. 

When cr^o-p takes on its maximum value for this model, the light Higgs becomes quite 
light lying just above the LEP 200 limit. Thus for tan/3=7 one finds 100 GeV< < 107 
GeV which would make the h relatively easy to find at the LHC and the Tevatron RUN 
II. In contrast the squark can become quite heavy. This can be seen from Eq. (9) where 
the negative nonuniversal term dominates the m\ contribution, and hence making Wq large 
reduces /i 2 , making cr^o_ p larger. Thus the first two generation squarks for this situation 
have masses in the range 600 GeV< nig < 1200 GeV, and lie above the gluino. Large values 
of rriq tend to suppress proton decay amplitudes of GUT models and would help relieve the 
tension between a large cr^o_ p and a small p decay rate. 



IV. D-BRANE MODELS 

Recent advances in string theory based on Dp-branes (manifolds of p+1 dimensions) has 
led to a revival of phenomenologically motivated string models. One class of such models 
making use of type IIB orientifolds |32| allows one to put the Standard Model gauge group 



on 5-branes, manifolds of six dimensions, of which four are usual Minkowski space and two 
are compactified on a torus. 

An interesting model of this class puts SU(3) C x U{1)y on one set of 5 branes, 5i and 
SU(2)l on a second intersecting set 52 [0. Strings starting on 52 and ending on 5i have 



massless modes carrying the joint quantum numbers of the two branes (presumably the SM 
quark and lepton doublets and the Higgs doublets), while strings begining and ending on 
5i have modes carrying SU(3)c X U(l)y quantum numbers (the right handed quark and 
lepton states). This leads to a model with a unique set of nonuniversal gaugino and squark 
and slepton masses at Mq parametrized as follows: 

rhi — rhz— —A = \/ r 3cos8i ) Qie~' tai m 3 /2 (10) 
m 2 = V3cos9 b (l - e 2 ) 1/2 m 3/2 



and 



m 12 



m 2 = (1 — 3sin 2 9b)m 2 



3/2sin 2 e b )m 2 3/2 (11) 

2 

.3/2 



Here rhi, i=l,2,3 are three gaugino masses at Mq , while m 2 2 are qi, II, Hi, H 2 sfermion 
and Higgs (mass) 2 and m 2 are Ur, d R , e R sfermion (mass) 2 . In addition, one has the B soft 
breaking mass and the /i parameter at the GUT scale 

B = \B \e i6oB ; fi = \fi o \e i0 ^ (12) 
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where ai, 9qb and 0$^ are possible CP violating phases. As can be seen from Eqs.(9,10), the 
parameters 0& and 0i are restricted to sinO b < l/v3 (9b < 0.615) and 0i < 1. 

Models of this type are of interest in that they are natural possibilities in string theory, 
and yet it would be difficult to see how such a symmetry breaking pattern could arise in 
conventional SUGRA GUT models. One aspect of the model is that the cancelations between 
the gluino and neutralino CP violating phases in the electron and neutron electric dipole 
moments (EDMs) implied by Eq. (10) allows the CP violating phases to be larger than usual 



at the elctroweak scale and still satisfy the experimental EDM bounds ||33|| . However, these 
experimental bounds combined with the requirement of radiative breaking of SU(2) x £7(1) 
at the elctroweak scale, leads to serious fine tuning of parameters at the GUT scale unless 
tan/3 ~ 3 — 5 |34| . Since we will be interested here in large tan/5 (to investigate the maximum 
values of o" x o_ p ) in the following we will first set all the CP violating phases to zero. One 
then has a model with four parameters and one sign i.e. 777.3/2, Ob, 6 and tan/5, and the sign 
of \i. (As before, fi must be negative for most of the parameter space to satisfy the b — > 57 
constraint.) 

Fig. H shows function of 777^,0, which is parametrized by 8b, ©1, 777 3 / 2 . One 

has that c^o_ p increases with increasing 6b since by Eq. (10) the squark masses of Fig. |2| 
decrease. Fig. |9] shows cr^o_ p for tan/3=15 and 777 3 / 2 =200 GeV. One sees that one needs a 

large tan/3, i.e. tan/3 ~ 15 to obtain a cross section within current experimental sensitivities, 
i.e. cr^o_ p ~ 1 x 10~ 6 pb. 

We next allow the CP violating phases to be non-zero to investigate their effect on 
cr^o_p cross section. We impose here the experimental constraint on the electron EDM of 
d e < 4.3 x 10~ 27 ecm at 95% C.L. [33]. (We do not impose the neutron EDM constraint 
here as there are a number of ambiguities in calculating d n p4]|). Fig. [U] shows that (T^o_ p 
decreases with increasing phase 2n — a>i, the presence of the phase decreasing & x o_ p , by 
a factor of two more. As mentioned above, such large phases, while still satisfying the 
experimental bound on d e , require a serious fine tuning of other parametrs. Fig. 11 shows 
the allowed range Ac/)^ of the phase <ft p to satisfy the electron EDM for tan/3=15, a.\ = 1.757T. 
One sees that M must be chosen very precisely to satisfy the EDM constraint on d e . If one 
were also to impose the neutron EDM constraint as well, the fine tuning would be even more 
severe [j35|. 



V. CONCLUSIONS 

We have considered here the question of what part of the SUSY parameter space can 
be probed by detectors sensitive to the neutralino proton cross section in the range of Eq. 
(1). In order to examine this we have considered the maximum theoretical cross section 
as function of 777^,0. The answer depends on the particular SUSY model one is considering. 

Thus for mSUGRA models, one requires tan/5 ~ 25 to achieve current detector sensitivities 
a xl-p > 1 x 10~ 6 pb. In this case, a large cr^o_ p corresponds to relatively heavy Higgs, 

rrih — (110 — 130) GeV, a light neutralino of 777^0 ~ 120 GeV (from the relic density 
constraint) and moderate squark masses (e.g. ~ (400 — 700) GeV lying below but close 
to the gluino mass. 

Nonuniversal SUGRA models can increase or decrease cr^o__ p by a factor of 10 to 100. In 
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the former case, one can begin to probe the parameter space for tan/3 ~ 4 for cr^o_ p > 1 x 10~ 6 
pb. Here the Higgs is relatively light, ~ (100 — 110) GeV, while the squarks are heavy 
(e.g. m rf - ~ 600 — 1200 GeV) and lie well above the gluino. Current data has in fact begun 
to eliminate part of the parameter space for tan/3 ~ 15. 

In the D-brane model considered, one requires tan/? ~ 15 for a^_ p > 1 x 10~ 6 pb. CP 
violating phases that can appear in such models lower the Xi ~~ P cross section by a factor 
of two or more. However, there is a serious fine tuning problem for the \i phase at Mq for 
such large values of tan /?. 

The fact that in the nonuniversal SUGRA models one can get a large cr%o_ p with a small 
tan/3 and large squark mass tends to relieve some of the tension between dark matter cross 
sections and proton decay. Thus if Xi dark matter were discovered with <J^o_ p ~ 1 x 10~ 6 
pb, such models could still have a low proton decay rate since r p is suppressed for small 
tan j3 and large rriq. 

In our analysis here, we have neglected the possibility of coannihilation effects since the 
rapid early universe annihilation produced by such effects generally require raising mo and 
mi/2 to avoid violating the lower bound of Eq. (6), lowering the Xi ~ P cross section . 
(We have been intersted here in the maximum values of a^o_ p .) We note, however, that 
aside from the coannihilation region arising from the f — Xi degeneracy for mSUGRA at 



low and intermediate tan/3 discussed in [ll|, there are three other coannihilation regions 
for the models considered here. In mSUGRA with large tan/?, the chargino xf can become 
degenerate with the x? leading to a new coannihilation domain. In the nonuniversal SUGRA 
model, one can see from Eq. (8) that the choice 5§ < can make the tr light and hence 
become degenerate with the Xi even for small tan/3. For the D-brane models one sees 
from Eq. (9) that as 0i gets large, m 2 and rh\ become equal allowing the Xi to become 
degenerate with Xi (roughly for Oi ~ 0.8). Unlike the case considered in [II], each of these 



new domains can occur for light neutralinos where cr^o_ p may be large. Coannihilation effect 



will be discussed elsewhere 36 
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FIGURES 




FIG. 1. Early universe annihilation of Xi through s channel Higgs (h, H, A) and Z poles and 
t channel squark and slepton (/) poles. 




FIG. 2. Scattering of x? by a quark of a nucleon through s channel q and t channel h, H, A, 
Z poles. 
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FIG. 5. mfj when cr^o_ p takes on its maximum value for tan/3 = 30 (top curve) and tan/3 = 50 
(lower curve). 
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FIG. 6. Maximum <7^o_ p for tan (3 = 7 for mSUGRA (lower curve) and nonuniversal model 
(upper curve). 
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FIG. 7. Maximum c^o_ p for nonuniversal models for tan/3 = 7 (upper curve) and tan/3 = 5 
(lower curve). 
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FIG. 8. for the D-brane model for m 3 / 2 = 175 GeV, tan/3=10 for 6^=0.5 (upper curve), 

#6=0.2 (lower curve). (The gaps are excluded regions.) 
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FIG. 9. a x°-p f° r ^ ne D-brane model for m 3 / 2 = 200 GeV, tan/3=15 for 9b=0.5 (upper curve), 
06=0.2 (lower curve). 




FIG. 10. a Xi~p ^ or ^b = ^-^' m 3/2 = 200 GeV, tan/3=15 for ct\ = 0, 1.85-7T, 1.757T in decreasing 
order. 
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FIG. 11. The allowed range A<^ of the phase 0^ for 0^=0.4, m 3 / 2 = 200 GeV, tan/3=15 for 

OL\ = 1.75-7T. 
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